Particle acceleration in relativistic current sheets 
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Relativistic current sheets have been proposed as the sites of dissipation in pulsar winds, jets in 
active galaxies and other Poynting flux dominated flows. It is shown that the steady versions of these 
structures differ from their nonrelativistic counterparts because they do not permit transformation 
to a de Hoffmann/Teller reference frame, in which the electric field vanishes. Instead, their generic 
form is that of a true neutral sheet: one in which the linking magnetic field component normal 
to the sheet is absent. Taken together with Alfven's limit on the total cross- field potential, this 
suggests plasma is ejected from the sheet in the cross-field direction rather than along it. The 
maximum energy to which such structures can accelerate particles is derived, and used to compute 
the maximum frequency of the subsequent synchrotron radiation. This can be substantially in 
excess of standard estimates. In the magnetically driven gamma-ray burst scenario, acceleration of 
electrons is possible to energies sufficient to enable photon-photon pair production after an inverse 
Compton scattering event. 
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Energy release by the process of magnetic reconnection 
has been suggested as the mechanism responsible for the 
production of energetic particles in the jets associated 
with active galactic nuclei (AGN) 0, 0, [1 II I! 13 and 
in pulsar winds H H ES El :t is thought that the 
same process may also be at work in gamma-ray bursts 
EJ 0, \uL uM- Reconnection has a large plasma 
physics, solar physics and magnetosphcric physics orien- 
tated literature, that is well-summarised in two recent 
monographs [171 Il8| . But the physical conditions in the 
high energy astrophysics applications mentioned above 
are significantly different, because relativistic effects are 
important. 

In this letter, it is shown that the generic form of a sta- 
tionary, relativistic current sheet is that of a true neutral 
sheet. Typical boundary conditions permit the sheet to 
be very extended along the magnetic field direction, but 
a simple argument limits the extent in the direction per- 
pendicular to both the sheet normal and magnetic field 
vectors. This finding is used to estimate the maximum 
possible energy to which particles can be accelerated by 
the DC electric field of the sheet. The results are ap- 
plied to pulsar winds, gamma-ray bursts and jets from 
active galactic nuclei, all of which may contain a mag- 
netically powered relativistic outflow. It is shown that 
the commonly used estimate of the maximum frequency 
of synchrotron radiation emitted by an accelerated elec- 
tron hv < 100 MeV can be substantially exceeded when 
Poynting flux dominates the energy flow. 

The current sheets at which reconnection and parti- 
cle acceleration takes place in astrophysics are relativis- 
tic in two senses: Firstly, the magnetization parame- 
ter, a (defined in Eq. QJ, is large and the Alfven speed 
va = cy/a/(l + a) is close to c. Secondly, the ge- 



ometry of the current sheet at which magnetic energy 
is dissipated is dictated by a highly relativistic plasma 
flow. Particle acceleration depends crucially on both 
the magnetization parameter and the field configuration. 
However, most analytic treatments utilize nonrelativis- 
tic models of current sheets. Numerical simulations of 
current sheets with a rs 1 have been performed using 
a two-dimensional PIC code 0], but the field geometry 
was constrained to be close to the standard Sweet-Parker 
configuration. Analytic work on relativistic reconnection 
[ill Efl Ell E^ ] has not yet addressed the question of 
particle acceleration. 

The relativistic effects associated with a large mag- 
netization parameter are readily appreciated. On the 
other hand, the geometrical effects of a relativistic flow 
are more subtle. The situation is closely analogous to 
that of MHD shock fronts, which can be classified into 
"subluminal" and "superluminal" according to whether 
the speed of the intersection point of the magnetic field 
and the shock front is less or greater than c 0] • In 
each Lorentz transformation enables the shock 

to be viewed from a reference frame in which it has a 
particularly simple configuration. 

Subluminal shocks permit a transformation to a 
de Hoffmann/Teller frame, where the electric field van- 
ishes both upstream and downstream. Nonrelativistic 
shocks, at which the speed /3 up of the shock front in 
units of c observed in the upstream medium is small, 
/3 U p <C 1, are subluminal, provided the angle between 
the magnetic field direction and the shock normal in this 
frame 6> up satisfies cos0 up > (3 up . Fine-tuning of the up- 
stream magnetic field direction would be required to vi- 
olate this condition and, given that no MHD wave mode 
can race ahead of the shock front, this seems unlikely. 



2 



FIG. 1: (Color online) The striped pattern of a pulsar wind. 
A magnetic dipole embedded in the star at an oblique angle 
to the rotation axis introduces field lines of both polarities 
into the equatorial plane. The current sheet separating these 
regions is shown. In the inset, an almost planar portion of 
this sheet (dashed line) is shown, together with the magnetic 
field lines, assuming they undergo reconnection. 
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On the other hand, the simplest configuration for a su- 
pcrluminal shock is one in which the magnetic field lines 
are perpendicular to the shock normal both upstream 
and downstream — an exactly "perpendicular" shock. 
This is the usual configuration for shocks which are rel- 

— 1/2 

ativistic, i.e., those with P = (l — /3jLj 3> 1. Then, 
fine-tuning of the upstream magnetic field direction such 
that sin 9 < 1/T would be needed to avoid superluminal 
speed of the intersection point. 

Thus, the generic nonrelativistic shock can be pic- 
tured as a stationary surface crossed by magnetic field 
lines along which the plasma flows. The generic rela- 
tivistic shock, however, can be seen as a stationary sur- 
face through which a magnetic field orientated precisely 
perpendicular to the surface normal is advected by the 
plasma. 

To understand the analogy with current sheets, it is 
best to consider a specific example: the relativistic MHD 
wind driven by a pulsar — a rotating neutron star in 
which a magnetic dipole is embedded with axis oblique 
to the rotation axis. Approximating the dipole by a 
split-monopole, an asymptotic solution is available, valid 
at large distance from the star J2||: the wind pulls out 
the field lines into a striped pattern j2(| , of wavelength 
27r/3rL, (t~l = c/fi is the light-cylinder radius, VI the rota- 
tion speed of the star) that propagates radially at a speed 
(3c close to c. Figure shows the equatorial plane where 
a thin helical current sheet separates regions of magnetic 
field that emerge from opposite magnetic poles. This is 
the most widely discussed scenario for a pulsar wind, al- 
though no solution exists connecting it to realistic bound- 
ary conditions at the stellar surface and alternative pic- 
tures involving currents which do not collapse into thin 
sheets have also been proposed |2^. In ideal MHD, the 
current sheet is a neutral sheet in which the magnetic 



field vanishes. The inset in Fig. ^ shows that field line 
reconnection implies a finite "linking" field component 
perpendicular to the sheet. The term "reconnection" can 
be misleading here, since the field line configuration is 
stationary in the corotating frame. However, plasma cer- 
tainly leaves these field lines, so that the process qualifies 
as reconnection under a more general definition [la, |28| . 

At a large distance r 3> rt, from the star, the sheet can 
be considered as locally plane. Then the linking compo- 
nent is in the z direction, and the reversing component 
in the ±x direction. The idealised steady reconnection 
picture discussed here implies an electric field E v in the y 
direction that causes the magnetized plasma to drift into 
the current sheet from both positive and negative z. A 
shearing component of the magnetic field in the y direc- 
tion can also be added, but this vanishes in the equator if 
the magnetic and rotational axes are orthogonal and will 
be ignored here, in which case the electric field remains 
along the y axis. 

In a stationary flow, reconnection balances the rate at 
which magnetic flux is drawn out from the star, imply- 
ing that the point at which a given magnetic field in- 
tersects the current sheet must corotate. However, far 
from the star, at r 3> tl, the corotation velocity is large 
compared to c and the relativistic configuration of Fig. ^ 
implies a drift speed in the sheet that is formally super- 
luminal: E y /B z 1. This same situation also arises in 
the case of extra-galactic jets, provided they have zero 
net poloidal flux i.e., they draw out magnetic flux from 
a central objectQ]. In contrast, the current sheets usu- 
ally considered in laboratory, solar and magnetosphcric 
applications arc assumed to be of finite length in the x 
direction and to have subluminal drift speeds at their 
edges. 

A two-dimensional, stationary configuration has E y 
constant and, for simplicity, it is usual to assume B z 
is constant too. Then it is possible to simplify the fields 
by performing a Lorentz boost: a subluminal sheet can 
be transformed by a boost in the x direction at speed 
cE y /B z into the "de Hoffmann/ Teller" frame in which 
the electric field vanishes. In this frame the orbits are 
complex, some being trapped in the sheet and others es- 
caping along the magnetic field lines after one or more 
passages through the plane z — 0. This is the standard 
configuration for the study of orbits in nonrelativistic cur- 
rent sheets 

mill. 

The particle energy is a constant of 
motion, since E = 0, and no acceleration results unless 
an additional scattering process is invoked. It is then easy 
to see that in any other reference frame the acceleration 
process can be regarded as one of reflection off the mag- 
netic field structure anchored in the sheet. This implies 
that ions are accelerated more efficiently than electrons, 
as observed in, for example, the geomagnetic tail [sij . 

Superluminal sheets, on the other hand, are trans- 
formed by a boost along x with speed cB z /E y into a 
configuration with vanishing linking field B z , but finite 
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whereas a cold electron-proton plasma gives 



Tmax ~ & 
7max « aM/m 



for protons 
for electrons, 



E y . (This is, in fact, a true neutral sheet: the field con- factor 7 max after acceleration is, following Alfven, 
figuration originally used by Speiser to discuss the 
orbits of nonrelativistic particles.) In this case all trajec- 
tories are trapped in the sheet and acceleration in the ±y 
direction ensues, according to the sign of the charge, as 
soon as a particle drifts into the acceleration zone around 
z = 0, where \E y /B x \ > 1. 

Particle acceleration in current sheets with finite B z 
has been extensively investigated 0, but the vanishing 
of B z has important implications. Most discussions of 
reconnection treat a Sweet-Parker or Petschek configu- 
ration in which the length of the current sheet in the 
x direction determines the dissipation rate. The linking 
field B z is then crucial for the determination of the spec- 
trum of accelerated particles, and, especially, the maxi- 
mum permitted energy 0, 0] , since it is responsible for 
ejection of particles from the acceleration zone (see also 
[ljj ). However, relativistic current sheets, such as that il- 
lustrated in Fig. can extend over large distances in the 
x direction, depending on the nature of the boundary 
conditions, and the linking field plays no role in ejecting 
particles. Instead, acceleration in a relativistic sheet is 
controlled by its finite extent in the direction parallel to 
the electric field E y . This is limited not by the boundary 
conditions, but by local parameter values, as described 
by Alfven j^l and Vasyliunas 

The maximum energy to which a particle can be ac- 
celerated by the DC field in a relativistic current sheet 
is determined by the product of the electric field E y and 
the maximum extent Ay in the y direction of the sheet. 
Alfven showed that E y x Ay must be finite by observing 
that charged particles which drift into the sheet through 
a surface of area Ax x Ay must leave it in opposite direc- 
tions along the y axis. The current I y in the y direction, 
part of which could be carried also by particles enter- 
ing the sheet along this axis, must therefore exceed the 
rate of inflow through Ax x Ay times the particle charge. 
Ampere's law, which relates the strength of the reversing 
component of the field to the current I y leads to the con- 
dition \qEy\Ay < Bl/{Aim±), where n± is the number 
density of particles of charge q outside the sheet. This 
argument, originally limited to stationary sheet configu- 
rations, can also be applied to evolving ones 

The magnetization parameter a can be defined in the 
local plasma rest frame as the ratio of the magnetic en- 
thalpy density to the particle enthalpy density (including 
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This quantity equals the ratio of the Poynting flux to the 
kinetic energy flux of a cold plasma observed in a ref- 
erence frame moving perpendicular to the magnetic field 
direction. Assuming the plasma consists of cold electrons 
and positrons, and that a 3> 1, the maximum Lorentz 



with M and m the proton and electron masses, respec- 
tively. It is interesting to note that in a plasma in which 
the magnetic field and particle rest mass are in rough 
equipartition (a w 1), the upper limit given by Eq. Q 
coincides with that quoted by Lesch and Birk [2|. How- 
ever, this situation arises only in relativistic plasmas. 
In the interstellar medium, for example, a ss 10~ 9 or 
smaller, in which case the upper limit on the energy gain 
reduces to Mv\. Standard estimates of the interstellar 
magnetic field and particle density (1 /iG, 1 proton/cm 3 ) 
imply that electrons can be accelerated, at most, to only 
mildly relativistic energies. In this case, and in solar sys- 
tem applications, direct acceleration by the DC field may 
be masked by particle acceleration in the turbulence fed 
by reconnection or the associated shocks [3(| . 

In addition to the limit imposed by finite Ay, radia- 
tive loss processes suffered during acceleration may also 
affect the maximum particle (and especially electron) en- 
ergy. In the absence of a linking magnetic field compo- 
nent, the particle trajectories, as first found in the non- 
relativistic limit by Speiser [32|, are attracted to the ac- 
celeration zone around the 2 = plane, where they un- 
dergo essentially one-dimensional acceleration in the elec- 
tric field. In this case the losses by synchrotron radiation 
(or, more precisely, radiation caused by the electromag- 
netic fields of the sheet) are unimportant. For large a, 
this enables synchrotron photons of relatively large en- 
ergy hv m&x w (heB /mc) Jmax ^o be generated when the 
accelerated particles leave the current sheet after travers- 
ing the distance Ay. This can substantially exceed the 
frequently quoted (and magnetic field independent) up- 
per limit of /ifmax ~ 100 MeV obtained by setting the 
synchrotron cooling rate equal to the inverse of the gyro 
period. For example, a relativistic current sheet in the 
wind of the Crab pulsar just outside the light cylinder 
(B « 10 6 G, a « 10 4 H2)' is capable of accelerating elec- 
trons which subsequently emit synchrotron photons of 
energy up to 50 GeV, providing a plausible production 
mechanism || for the flux of GeV gamma-rays observed 
by the EGRET experiment 37]. 

Unlike synchrotron radiation, inverse Compton scat- 
tering on an ambient radiation field is not suppressed 
for linearly accelerated particles, and can limit the max- 
imum energy, especially near a luminous source such as 
an AGN. Parameterising the accelerating electric field by 
V) = \E y /B x \ < 1, this limit can be written for a current 
sheet embedded in a jet that emerges from the AGN with 
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bulk Lorentz factor T = Tio/10 as 



In 



10 5 ( V L P47 ) 1/4 (r 10 R 12 /L hh47 ) 1/2 . 



(5) 



where Lp47 and I/bb47 are the luminosity carried by 
Poynting flux in the jet and the luminosity of the central 
photon source, respectively, in units of 10 47 ergs~ 1 , and 
R = 10 12 i?i2 cm is the distance of the sheet from this 
source. Equation (JHJ assumes that the Compton scatter- 
ing takes place in the Thomson limit, which is valid for 

/max rsj lu L • 

In the physical conditions thought prevalent at the 
photosphere of a magnetically accelerated gamma-ray 
burst outflow [3 {Lp47 « 10 3 , L hhi7 w 50 and T w 100 
at the photosphere i? = 10 11 cm), Eq. J5J indicates 
a maximum electron Lorentz factor (in the comoving 
frame) of 7 max ~ 4 x 10 5 . However, since the photon field 
in this case has a relatively high temperature (~ IkeV), 
we conclude that Compton scattering is unable to limit 
the maximum electron energy to a value such that the 
Thomson limit applies, and, consequently, that photon- 
photon pair creation is likely to ensue. A feedback mech- 
anism appears possible here, since pair creation will de- 
press the magnetization parameter er, and, hence, the 
limit given by Eq. J3J|. 

Relativistic current sheets are known to be unstable 
to the growth of the tearing mode and other insta- 
bilities are also likely to operate (see, for example, [39|). 
On scale lengths comparable to the sheet thickness it is 
likely that an oscillating component of B z will be gener- 
ated. As discussed in standard nonrelativistic pictures, 
this may, on small scales, cause breakup of the sheet into 
a series of magnetic islands and X-points arranged alter- 
nately along the x-axis, providing stochastic deflections 
of the particle trajectories within the sheet. Such a scat- 
tering mechanism is, in fact, essential for the formation 
of an equilibrium structure such as the Harris sheet |40| , 
and it may also play a role in enhancing the effect of syn- 
chrotron losses. However, in a steady relativistic sheet 
the boundary conditions do not permit this scattering 
process to evacuate the plasma from the dissipation re- 
gion, because they prevent the formation of a stationary 
pattern with outflow along the a; axis. On the other hand, 
a quasi-stationary pattern containing relativistic sheets 
of finite extent in the y direction does appear possible, 
and enables estimates of the maximum permitted energy 
of particles accelerated by the DC field to be made. In 
reality, relativistic sheets may dissipate in a highly non- 
stationary fashion as suggested in the nonrelativistic case 
by, for example, solar observations |36| . but clarification 
of this must await the results of numerical simulations of 
the generic relativistic sheet configuration. 

I thank Y. Lyubarsky for helpful comments on the 
manuscript. 
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